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f~] ' Abstract. The interstellar medium heated by supernova explosions (SN) may acquire an expansion velocity larger 

^ , than the escape velocity and leave the galaxy through a supersonic wind. Galactic winds are effectively observed 

• in many local starburst galaxies. SN ejecta are transported out of the galaxies by such winds which thus affect 

' the chemical evolution of the galaxies. The effectiveness of the processes mentioned above depends on the heating 

, efficiency (HE) of the SNe, i.e. on the fraction of SN energy which is not radiated away. The value of HE, in 

' particular in starburst (SB) galaxies, is a matter of debate. We have constructed a simple semi-analytic model, 

^ I considering the essential ingredients of a SB environment which is able to qualitatively trace the thermalisation 
history of the ISM in a SB region and determine the HE evolution. Our study has been also accompanied by 
fully 3-D radiative cooling, hydrodynamical simulations of SNR-SNR and SNR-clouds interactions. We find that, 
as long as the typical time scale of mass-loss of the clouds to the ambient medium, which is often dominated by 

I ■ photoevaporation, remains shorter than the time scale at which the SNRs interact to form a superbubble, the SN 

' heating efficiency remains very small, as radiative cooling of the gas dominates. If there is a continuous production 

, of clouds by the gas swept by the SNR shells, this occurs during the first <16 Myrs of the SB activity (of ~ 30 

' Myrs), after which the efficiency rapidly increases to one, leading to a possible galactic wind formation. Under an 

I extreme condition in which no clouds are allowed to form, other than those that were already initially present in 

II , the SB environment, then in this case HE increases to one in only few Myrs. We conclude that the HE value has 

Q ' a time-dependent trend that is sensitive to the initial conditions of the system and cannot be simply assumed to 

^ , be ~ 1, as it is commonly done in most SB galactic wind models. 
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1. Introduction In a SB region, several SN explosions occur at a high rate 

1 / n \ 1 inside a relatively small volume, and their remnants are 

Energization of the interstellar medium (ISM) by super- n i ^ • j. j. -j-i, u j-i, t tTn'7A\ • i,- 

° . , , . , likely to mteract with each other. Larson (1974), m his 

nova (SN) explosions covers a crucial role in a number . . i rxi • i r ±- 

^ ' ^ pioneering work; oi the monolithic galaxy lormation, com- 

of astrophysical situations. During the galaxy formation, ^^^^^ ^^^^^^^ ^^^^^^ ^^^^^^ ^ 

the first exploding stars may heat and expand the lo- , rn ^ r i i u j- j. j 

, ,11 r when 50 per cent of its thermal energy has been radiated 

cal gas reducing, or even halting the ongoing star forma- , cato • j. j.- j.- j- j c j j-i, 

. . „ . . away, to the SNRs interaction time, tint, defaned as the 

tion. Galaxies with highly active star formation regions, . i r ii j. j. nn x. 

° 1 , , time it takes for the remnants to occupy ~ 60 per cent 

starburst (SB) galaxies, present a high SN rate. The gas p , £ j.u ct> • Ks:^ ■ i ^- ^■ 

, . „ ,01 the volume of the SB region. After an interaction time, 

heated by SNe may acquire an expansion velocity larger , , n-j j u j j-u ■ 

J H f J b |-j^g remnants have collided with each other and their evo- 

than the escape velocity and leave the galaxy through a , , • d i- ; • i • --f + ^ + 

, „ , , 1 , lution ceases. Radiative losses are important if tc < tint- 

supersonic wind. Galactic winds are enectively observed j i ^ j-i, ^ • 4.1, c ■ x. j. j.v, 

^ 1 /I Nr. Larson evaluates that, m the cases of interest, the energy 

m many local SB galaxies (Lehnert & Heckman 1996). SN 1 • 1 • a- 1 ^ -ij. j j. ii, • m ^ 

, r , , , which IS eftectively transmitted to the gas is ~ 10 per cent 

ejecta are transported out of the galaxies by such winds r , . nir,- 1 £4.1, ttt:^ u u 

rr- , ■ , ? • r 1 i • explosiou eucrgy. ihis value of the ML has been 

which thus aftect the chemical evolution of the galaxies. j ^ j • u r cm i 1 c n 1 c n 

° adopted m a number of SB models (e.g. Babul & Rees 

The effectiveness of the processes mentioned above de- ^993, Murakami & Babul 1999). Bradamante, Matteucci 

pends on the heating efficiency (HE) of the SNe, i.e. on ^ D'Ercole (1998), following a different line to estimate 
the fraction of the SN explosion energy that remains ef- ^^^^ considered only the final fraction of the SN en- 

fectively stored in the I SM gas and is not radiated away. ^^^^ ^^^^ effectively turns into thermal energy of the am- 

Send offprint requests to: Melioli C. bient medium when the remnant, already at the end of its 
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radiative phase, stalls. This gives HE ~ 0.03 in the cases 
of interest, and this value has been adopted by Ferrara 

& Tolstoy (2000) and by Rccchi, Mattcucci and D'Ercole 
(2001) in their models of dwarf galaxies. 

These low values of the HE have been questioned by 
several authors. Following Cox & Smith (1974), Larson 
(1974) warned that, for a high SN rate, the collisions of 
the first generation of SNRs possibly lead to the forma- 
tion of an interconnected network of tunnels containing 
very hot, tenuous gas maintained at high temperature by 
continuing supernova explosions. The successive genera- 
tions of SNRs are likely to produce a value of the HE close 
to unity. In fact, the remnants do not reach high densi- 
ties during their expansion, and their radiative losses are 
expected to remain negligible as the emissivity is propor- 
tional to the square of the density. Dekel & Silk (1986) 
find that the SNRs interact before entering in the ra- 
diative phase, and almost all their energy is transferred 
to the ISM gas. Short interacting times seem to invali- 
date the assumption of stalling remnants by Bradamante, 
Matteucci & D'Ercole (1998). Strickland & Stevens (2000) 
argue that with a value of the HE as low as that proposed 
by Bradamante, Matteucci & D'Ercole (1998) the galactic 
wind observed in M82 (the best studied SB with a galactic 
wind) could be hardly driven. Chevalier & Clegg (1985) 
find that 100 per cent is a likely value for the HE in M82. 

In the absence of a clear response from observations 
(e.g. Delia Ceca et al. 1999; Cappi et al. 1999), the above 
considerations have led to the widespread belief that the 
HE within a SB region must be quite high to sustain galac- 
tic winds. For this reason, most of the simulations of galac- 
tic winds found in literature assume a HE value of 100 per 
cent (e.g. Suchkov et al. 1994, Silich & Tcnorio-Tagle 1998, 
D'Ercole & Brighenti 1999, Mac Low & Ferrara 1999, 
de Gouveia Dal Pino & Medina Tanco 1999, Strickland 
& Stevens 2000). However, several of these "standard" 
models fail in reproducing the sizes of the active regions, 
which turn out to be larger than those observed (see, 
e.g. Tcnorio-Tagle & Munoz-Tuhon 1998 and references 
therein). Moreover, Recchi, Matteucci & D'Ercole (2001) 
(also Tosi 2003), are able to reproduce the chemical and 
dynamical characteristics of the SB IZwlB, for example, 
only assuming, in an instantaneous SB, a HE value of 0.03, 
which is about 30 times lower than the value commonly 
assumed in the simulations of SB galactic winds. In mod- 
els with HE=1, all the gas in IZwl8 is lost, in opposition 
to what is observed. We are thus motivated to investigate 
further the coupling between SNRs and the ISM inside ac- 
tive regions, and explore under which circumstances the 
HE may assume low values. In order to do this, we will not 
be concerned with a detailed description of the structure 
of the star-bursting galactic region. Instead, we construct 
a model, as simple as possible, disregarding all the details 
that are not essential for our purposes; such a model is 
highly idealized, but still able to give us insights on the 
thermalisation history of the ISM inside a star-bursting 
region. 



In the following sections, we outline the characteris- 
tics of the model and present solutions for a steady state 
system (§ 2); discuss the physical processes in a SB envi- 
ronment including the presentation of the relevant results 
of 3-D numerical simulations of the interactions between 
SNRs and clouds (§ 3); and study the relevant time scales 
for the mixing of the several gaseous components of the 
system (§ 4). In § 5, we present the results of our evolution 
model for a SB environment including all the physical pro- 
cesses described in the previous sections; and finally in § 
6, we present a brief discussion and draw our conclusions. 

2. The model 

2.1. Assumptions 

The events occurring in a SB region are quite complex. 
Supernova remnants (SNRs) expand through a diffuse 
medium in which clouds are embedded. These clouds suf- 
fer mass loss by ablation, thermal evaporation and photo- 
evaporation due to the huge amount of ionizing photons 
emitted by massive stars and by SNs. The gas is heated by 
such SNs, and is cooled via radiative losses. These losses 
depend on the square of the gas density, so that two com- 
peting facts occur: on one hand the SNRs sweep the diffuse 
gas, thus lowering its mean density in the SB region; on 
the other hand the gas density is increased by mass loss 
from the clouds. The importance of the radiative cool- 
ing depends on which of these two processes dominate. If 
the cloud mass loss rate is smaller than the sweeping rate 
by SNRs, the successive SNRs expand in a medium more 
and more rarefied, thus remaining adiabatic and eventu- 
ally transferring most of their energy to the ISM; in this 
case the HE would be close to unity. On the contrary, HE 
would be very low for a cloud mass loss rate larger than 
the sweeping rate. 

It is very difficult to make numerical simulations of 
the above scenario taking into account all the physical 
processes. Moreover, the length scales involved are quite 
different, ranging from hundredths of parsec of the con- 
duction front width at the edge of the clouds (see next 
section) to the size of the SB region, typically 100 pc 
(e.g. Meurer et al. 1995, Planesas, Colina & Perez-Olea 
1997); numerical simulations able to cover this range of 
length scales with appropriate spatial resolution would 
require an unrealistic amount of cpu time. We thus fol- 
low a semi-analytical approach, as simple as possible, but 
incorporating the "micro" physics as found in literature. 
We refer to previous works that have addressed the evo- 
lution of a multiphase ISM assuming realistic cooling and 
heating processes and a number of mechanisms for clouds 
destruction and formation (see, e.g., Bertoldi & McKee 
1990, Rosen & Bregman 1995; Vazquez-Semadeni, Pasot 
& Pouquet 1995; Shore & Ferrini 1995; Wada & Norman 
1999, Wada & Norman 2001), though they were mostly 
concerned with the ISM of our Galaxy and have focused 
on regions with typical sizes and physical conditions which 
are generally distinct from those addressed here, most of 
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these studies have provided the guide hues for building 
the present approach. 

Wc consider the following toy model. An instantaneous 
burst of star formation occurs inside a spherical region of 
radius i?sB- Assuming the Salpeter initial mass function 
(IMF) the number M of stars with mass greater than SAf© 
(i.e., the SN number) is A/" ~ 0.01(Mb/MQ), where Mb is 
the total mass of the stars in the SB. The type II supernova 
(SNII) activity lasts up to a time % ^ i x 10^ yr, which 
is the lifetime of an 8 Mq star. A constant SN rate is 
thus given by 7?. = N /%, which is in good agreement with 
more accurate evaluations given by Leitherer ct al. (1999). 
It follows that the mass deposition rate and the energy 
injection rate are given, respectively, by Msn = 107^ Mq 
yT~^ and _Esn = \O^^Tl. erg yr~^, where wc have assumed 
an energy output of 10^^ erg and a mass return of lOM© 
per SNII. 

Assuming a star formation efficiency of ^ 10 % 
(Colina, Lipari & Macchetto 1991) and similar character- 
istics to those observed in the Rosette cloud, where 78 % of 
the gas is in the form of compact globules (Williams, Blitz 
& Stark 1995), we obtain that the total mass of clouds in 
the SB region is Mg = 8Mb. 

The diffuse gas is assumed to be optically thin, and 
the cooling function is approximately A(T) = 1.6 X 
10-i9/3T-°-5 erg cm^ s'^, for 10^ K < T < 10^-^ K 
(McKee & Begelman 1990), where the factor /? is one if the 
gas is in ionization equilibrium. At the conductive inter- 
faces, non-equilibrium ionization effects can increase the 
cooling up to an effective value /3 ~ 10 (Borkowski. Balbus 
& Fristrom 1990, McKee & Ostriker 1977). Following 
Mathews & Brcgman (1978), in the range 10* K < T < 
10^ K, wc assume a linear dependence for A(T) with T, 
A(T) = 5.35 X 10"27y ^j.g pj^3 g-i 

As a further simplification, we assume that all the 
quantities are uniformly distributed inside the spherical 
SB volume, and that the gas leaves this region with a 
Mach number Ai = 1, i.e. the gas moves with a veloc- 
ity equal to its sound speed (cf. Chevalier & Clegg 1985), 
corresponding to a free expanding steady state wind ex- 
tending to infinity. A schematic picture of the model is 
shown in Figure 1 and Table 1 summarizes the typical 
values of the parameters of a SB environment. 

The aim of the model is to calculate the HE defined as 
the ratio of the enthalpy flux through the boundary of the 
SB region to the SN luminosity (see cq. (29)). As a first 
step, we present in the next subsection possible steady 
solutions under very simplistic assumptions. Despite its 
excessive idealization, this study outlines the role of the 
mass return from the clouds. 

2.2. Steady state solutions 

In this simple example, we do not deal with the physics 
of SNR expansions and of mass loss from the clouds, but 
simply consider SNs and clouds as sources of energy and 
mass. Assuming that the gas moves with a velocity equal 




Fig. 1. Schematic picture of the model: a spherical SB 
with radius Rsb hosting expanding SNRs and clumps of 
dense and cold clouds 

to its sound speed, mass and energy conservation then 
yield: 

dp 3 . , 

where p and p arc the density and the pressure of the 
diffuse ambient gas, psN represents the mass injection rate 
by SNe, pa is the mass loss rate by the clouds due to their 
destruction, mu is the proton mass, p = 1.3 is the mean 
mass per nuclei of the ionized gas assumed with 90 % H 
and 10 % He abundances, = (2.1fc//imH)^/^TV2 is the 
isothermal sound speed, Rsb is the radius in parsec of 
the SB, esN is the rate of energy injection by the SNe 
explosions. Vc is the mean cloud velocity; at this stage 
we do not make any assumption about the mechanisms 
responsible for the cloud mass-loss, and simply assume 
Vc — 10^ cm s~^ as a fiducial value. 

In order to obtain possible steady solutions of the equa- 
tions (l)-(2), we drop the time derivatives and, after some 
simple algebra, we obtain 

3(1 + Ocl + %Psn(1 + enc.) = + (3) 
where ^ = Pci/psN, k is the Boltzman constant, and 

T{ \-i^Cs~^ ifr<105K 
^""'''XKc-^ ifT>105K 

Here K = 1.6 x 10-1^(2. lfc)i/2(^TOH)-^/^ cm^ s"* g-^ and 
Cg is the sound speed at T = 10^ K. 

The l.h.s. of equation (3) describes the gas cooling. In 
particular, the first term represents the cooling due to the 
adiabatic expansion, while the second term represents the 
radiative losses. The r.h.s. contains the heating sources. 
The second term represents the SN heating per unit mass 
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Table 1. Values for the characteristic parameters of a SB galaxy assumed in this study 



Model I RsB (pc) I Mb (Mo) I Mg (Mo) I | 7^ (yr-^) | h (Myr) 

Spherical I 100-700 I 10^ - lO'' I 10® - 10^ I 10^ - 10^ I 3 x 10"^ - 3 x 10"* I 30 
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Fig. 2. Solid lines: cooling curves for /3 = 1 and three dif- 
ferent values of ^. Dashed lines: the same for /3 = 10. The 
heating term is represented by the horizontal line which 
does not depend on /3 and is virtually independent of ^. 

of the ejecta. We note that this term is an intrinsic prop- 
erty of the SNe and is independent of any assumption 
about the SB. The first term represents the thermalisa- 
tion of the kinetic energy of the gas lost by the clouds 
moving with a mean velocity v^. We note that this term 
remains negligible with respect to the second one for any 
reasonable value of ^ and v^. 

Figure 2 illustrates the heating sources and the energy 
sinks of equation (3) as functions of T for a model with 
i?SB = 100 pc and Mb = 10^ Mq. The three solid curves 
represent the l.h.s. of equation (3) for three different val- 
ues of ^ and for /? = 1; the dashed curves hold for the same 
values of ^, but for /3 = 10. The horizontal line represents 
the r.h.s. of equation (3) and is virtually independent of ^. 
At high temperatures, radiative losses are negligible and 
the gas cools via adiabatic expansion (the rising branch of 
the cooling curve in Fig. 2). The opposite is true at inter- 
mediate temperatures, and the radiative cooling increases 
with decreasing temperature. For T < 10^ K, cooling is 
constant and independent of the temperature. 

In Fig. 2 ^ , three different cases are recovered, depend- 
ing on the value of ^. For ^ < (see below) only one 

^ We notice that due to the uncertainties associated to the 
value of /3 in the radiative cooling curves of Fig. 2 (e.g., Cowie, 
McKee, & Ostriker 1981, Borkowski, Balbus & Priston 1990), 



steady solution exists (the intersection point between the 
heating line and the cooling curve in Fig. 2) and it occurs 
at a high temperature (the hot solution). As pointed out 
above, at high temperatures radiative losses are negligible, 
and the HE value is close to 1. Following Field (1975), it is 
easy to show that this solution is stable. For > C > ^h, 
two sohitions are possible, the hot solution and the warm 
solution. This latter solution has very low values of HE 
but is highly thermally unstable. Numerical integration 
of equations (l)-(2) shows that the system moves toward 
higher or lower temperatures depending on the value of ^ 
(of course, lower values of ^, and hence of density, favour a 
decrease of the cmissivity and an increase of the temper- 
ature). Finally, for ^ > no steady solution is possible, 
and the system moves toward low values of temperature 
and HE, for any reasonable choice of initial conditions. 
For the model shown in Fig. 2, we have = 17.5 and 
= 70.5. For this model Tl = 10"^ yr and the rate of 
mass deposited by SNe is Msn = 10"^ Mq yr"^ Thus, if 
the rate of mass transfer from the clouds to the diluted gas 
is larger than ~ 0.7 Mq yr^^, the system cannot seat on a 
steady solution and the HE is expected to be very low. For 
cloud mass transfer rates smaller than Msn = 10~^ Mq 
yr~^ but larger than ~ 0.17 Mq yr~^, both evolutions 
toward high values or toward low values of the HE are 
possible. It is noteworthy that, for a SN rate 100 times 
larger, as that estimated for M82, we have Msn = 1 Mq 
yr~^ and ^„ = 14.3. In order to obtain low values of the 
HE, clouds should lose mass at a rate greater than 14 Mq 
yr~^. As shown below, such a high rate can be sustained 
only for short times compared to the SB age, ti,. The value 
of the HE is thus expected to be close to unity in this case, 
as suggested by Chevalier & Clegg (1985). 

As shown by this simple example, the behavior of the 
ISM in the SB region is regulated by the value of ^. Thus, 
in order to make a more realistic model, we must take into 
account in more detail the mass exchange mechanisms be- 
tween clouds and diffuse ISM. They include SNRs inter- 
actions, clouds-SNR interactions, photoevaporation, drag 
due to Kelvin-Helmholtz instability and thermal evapora- 
tion. We discuss these mechanisms in the next section. 

3. Physical processes in the SB environment 

3.1. Clouds formation 

Cold dense gas is observed in SB regions, as in M82 (Kerp, 
Herbstmeier & Mebold 1993; Cecil et al. 2001; Pietsch 
et al. 2001 ). Radio observations of SB galaxies indicate 

we have adopted /3 ~ 1 in all the calculations hereafter, as it 
gives a lower limit condition for the gas cooling (see below). 
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the presence of clouds with masses between lO^M© and 

5 X IO^Mq, radii between 0.5 pc and 1 pc, and tempera- 
tures between 50 K and 200 K (e.g. Ccsaroni ct al. 1991; 
Carral et al. 1994; Paglione, Tosaki & Jackson 1995; Garay 

6 Lizano 1999). However, we could also expect the pres- 
ence of even smaller clouds which arc not detectable due 
to observational limitation. In our own Galaxy, molecular 
clouds are detected in a range between 1 Mq and lO^M© 
with a mass spectral index distribution a = 1.5 (Scalo & 
Lazarian 1996). These molecular clouds show in turn sub- 
structures called clumps, that are distributed in a mass 
range between 10~^ and 10^ M©, with power law spectral 
index for the distribution in a range a = 1.4-1.9 (Kramer 
et al. 1998; Blitz 1993). 

While most of the cold gas observed in SB regions is 
pristine, new cold clouds may form under the action of 
SNRs. CiofR and Shull (1991) have shown that many small 
clouds (with sizes of fractions of parsec) form after the in- 
teraction among SNRs generated by a high rate of SNe 
explosions. Scalo & Chappell (1999) highlight the forma- 
tion of cold filaments in HIT regions, as a consequence of 
strong stellar winds and SNe explosions occurring in a re- 
gion containing a large number of stellar sources. 

Wc have recently run three-dimensional (3D) hydro- 
dynamical simulations of SNRs interactions employing a 
modified version of the Yguazu-a adaptive grid code orig- 
inally developed by Raga, Navarro-Gonzalez & Villagran- 
Muniz (2000; see also Raga et al. 2002, Masciadri et al. 
2002, Gonzalez et al. 2004). This code integrates the hy- 
drodynamic equations explicitly accounting for the ra- 
diative cooling together with a set of continuity equa- 
tions for several atomic/ionic species employing the flux- 
vector splitting algorithm of Van Leer (1982). The follow- 
ing species have been considered: H I, H II, He I, He II, He 
III, C II, C III, C IV, O I, O II, and O III. The reaction 
rates and the non-equilibrium cooling function are given 
in Raga et al. (2002). The calculations were performed on 
a four-level, binary adaptive grid with a maximum resolu- 
tion along the x, y, and z axes of 0.1 pc. The computational 
domain extends over (25) x (3.125) x 25 pc, corresponding 
to x256 X 32 X 256 grid points at the highest resolution 
grid level (see also, Melioli, de Gouveia Dal Pino & Raga 
2004 ). 

The simulations have shown that the formation of very 
dense filaments (1000 times denser than the hot gas) and 
dense surfaces (100 times denser than the hot gas) is pos- 
sible in times <10^ yrs (Figure 3). Rayleigh- Taylor (R- 
T) and Kelvin-Helmholtz (K-H) instabilities arc able to 
fragment these structures and produce new clouds, in a 
process that can continue during all the SN activity. ^ 

^ An alternative mcehanism for clouds formation could be 
due to ISM turbulence (e.g., Ehlerova et al. 1997, Palous, 
Wunsch & Ehlerova 2001, Ortega, Volkov & Monte-Lima 2001, 
Elmegreen 2003). Dynamical instabilities and stellar pressures 
can drive turbulence in the ambient medium that in turn may 
cause the formation of clouds and filaments. However, as most 
of these studies relying on turbulence have neglected the pres- 
ence of SN shell interactions at a high rate, the resulting 
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Fig. 3. Grey-scale map of the density distribution showing 
the interaction of 5 SNRs in a box of 25 pc x 3.125 pc x 
25 pc, corresponding to 256 x 32 x 256 grid points at the 
highest grid level. The left and bottom faces that define 
the beginning of the computational box have reflective 
boundaries and the right and top faces have continuous (or 
outflow) boundaries. The ambient medium has n=l cm~"^ 
and T= 10^ K. The initial energies of the SN explosions 
are assumed to be; Eq = 6.3 x 10"*^ erg, 1.6 x 10"''" erg, 
4.2 X lO'^^ erg, 3.8 x lO^'^erg, and 1.3 x 10^° erg. The panels 
correspond to two different snapshots of the evolution of 
the system: t ~ 2.5 x 10^ yr (upper panel) and 4.5 x 10'' 
yrs (lower panel) (see also Melioli, de Gouveia Dal Pino 
& Raga 2004, for details). 



We assume that the rate of clouds formation depends 
on the instantaneous amount of total mass of gas in the 
SB region. We can express the rate of formation of the 
new clouds as: 

dt rrictint 

timescales for clouds formation and their resulting sizes are 
larger than those expected in the SB environment. 
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where Nc^int is the total number of clouds formed by SNRs 
shell interaction, Mg^h is the total mass of gas that fills 
the SB volume, nic is the cloud mass, ti^t is the SNR 
interaction time (see also § 3.4 below), and 5 is the fraction 
of the gas that is deposited in a SNR shell. Models for 
the evolution of a SNR (see, e.g.. Chevalier 1974) indicate 
values for 5 between 0.9 and 0.99. 



3.2. Photoevaporation 

The SB volume is filled by a huge amount of ionizing 

photons which arc produced by the massive stars. In or- 
der to estimate the number of ionizing photons in our 
SB system, we have run the "mappings" code developed 
by Leitherer et al. (1999), which is available in the site 
www. stsci.edu/ science/ starhurst99/ mappings. Using ini- 
tial conditions which are appropriate to our SB environ- 
ment, we have found that the SN phase begins at a time of 
3.61 xlO^ yrs, that we have taken as the zero-time in our 
model. At this time, the number rate of ionizing photons 
]s S = 2.25 X 10^^s~i, and after this time it decreases as 

Initially, when a cloud becomes exposed to the ion- 
izing radiation, it implodes and injects part of the mass 
into the ambient medium (Bertoldi 1989). A compressed 
neutral globule remains, which continues to evaporate by 
photoionization (Bertoldi & McKce 1990). 

If we take into account the presence of a non-zero mag- 
netic field in the cloud, this can increase its effective pres- 
sure and make it more resistant against photoevaporation. 
In fact, in the case of the ISM of our own Galaxy, the pres- 
sure balance of the clouds with the hot high pressure, dif- 
fuse component is often dominated by the magnetic pres- 
sure rather than the thermal pressure of the cloud (e.g.. 
Cox 1995, Bower et al. 1995). In the case of a magnetically 
dominated cloud, it can be shown that when the pressure 
exerted by the ionizing front balances the pressure of the 
neutral globule cloud, then (Bertoldi & McKee 1990): 

1/21 1 



Finally, it can be shown that the radius of the cloud 



2.75 X 10" 



1000/ &2/7 

where m' is the cloud mass in Mq, ric is the initial cloud 
gas density and the parameter b gives the initial magnetic 
field expressed as Bq = hn'J'^ fi G, and = Sig/B^ss"^ 
with = 5/10*^ s~i. In our reference SB model (see 
below) we obtain = 0.225 pc~^ s~^. For a magnetically 
dominated cloud b > 0.3, and the photoevaporation rate 
can be written as a function of the cloud mass as (Bertoldi 
& McKee 1990): 

Mjjy = -1.79 X 10-^^^^/^m'*/^ Mq yr^^ 



(7) 



where is a dimensionless factor: 

2.7 66/7 



1/7 



where Uc is given in cm~ . In this work, in the presence 
of photoevaporation, we have Uc ^ 10^ — 6 x 10^ cm"^^ 
which implies tp = 0.5 — 1 (see §) for b^l. 



is: 



rc = 2.6 X 10-^:^-1/7^8/21 



(8) 



It is worth noting that the O and B stars are able to 
photoionize the clouds only if these are closer than the 
Stromgren radius 



i?,= 66.9(^)1/3 pc. 



(9) 



Assuming a number of O and B stars ~ A/'/3, the filling 
factor of the photoionizated SB volume, ffs, is: 



//, = 3 X 10^ 



3R% 



(10) 



SB 



Thus, the cloud mass- loss rate due to photoionization 

must be scaled by the factor //s, if ffs < 1. We will see 
below, however, that f fs increases to 1 in a quite short 
time (< 2 X 10^ Myr) compared to the duration of the SN 
activity. 

3.3. Drag 

The cloud motion through the medium leads to mass-loss 
by the clouds due to the Kelvin- Helmholtz (KH) instabil- 
ity (e.g. Lin & Murray 2000). Defining x = pd p as the 
gas density ratio between the cloud and the diffuse ISM, 
the growth timescale of the K-H instability is given by 
(Klein, McKee & Colella 1994) tkh i^cX^^V^c Thus, 
the mass-loss rate mc/rKH can be written as 



Md ~ 1.3 X IQ-'^nx^'^lv^fi Mq 



(11) 



where v^fi is the cloud velocity in units of 10 cm s- and 
Tc is expressed in pc. Taking into account equation (8), we 
may express this rate as a function of the cloud mass: 

Md ~ 8.8 X lQ-^°nx^'^v^,eJ^-^/^mf'^^ Mq yr-^. (12) 

where x is obtained assuming that the outer shell of the 
clouds is at a temperature ~ 10^ K. 

3.4. SNR-clouds interactions 

Besides the classical drag, clouds may also be destructed 
by interaction with SNRs. This is a rather complex pro- 
cess, and its investigation deserves numerical simulations. 
Works in literature dealing with this problem are rare. 
Poludnenko, Frank and Blackman (2002) have recently 
performed 2-D simulations of the interaction of a strong, 
planar shock wave with a set of dense inhomogeneties, em- 
bedded in a more rarefied ambient medium. Their results 
show that after compression, re-expansion and destruc- 
tion, a mixing phase occurs that causes a high increase of 
the ambient density. 

More recently, motivated by the investigation in the 
present work, we have carried out fully 3-D hydrodynam- 
ical radiative cooling simulations of the interaction of a 



Melioli et al.: SN heating efficiency of the ISM of SB 



7 



SNR shell with dense clouds embedded in a hot gas. As de- 
scribed before, the simulations were performed with an un- 
magnctizcd version of the adaptativc grid code YGUAZU 
(Raga, Navarro-Gonzalez & Villagran-Muniz 2000; Raga 
et al. 2002; Masciadri et al. 2002) which solves the gas dy- 
namical equations together with a set of continuity equa- 
tions for several atomic/ionic species without including 
the effects of photoevaporation and thermal evaporation of 
the clouds. With the radiation rates and non-equilibrium 
cooling function given by Raga et al. (2002) (see also 
Melioli, de Gouveia Dal Pino & Raga 2004 for details) 
our results, somewhat distinct from those of Poludnenko, 
Frank and Blackman (2002), show that after a time be- 
tween 10^ and 10^ yr the clouds are largely fragmented. 
However, little mixing occurs between the mass lost by the 
clouds and the ISM; instead, we find that new filaments 
form and the shell mass grows if the initial density of the 
clouds is of the order of that of the shell (see Fig. 4) . Only 
in few cases, when there are several very little clouds (with 
rc ^ hsh and ric > rish, where hsh and Ush are the thick- 
ness and the number density of the shell, respectively) the 
mixing occurs and the ISM density effectively grows by a 
factor 10-100 times larger than that of the hot gas density. 
Thus, as a rule, this type of interaction does not affect the 
diffuse ISM global density but contributes to form more 
smaller clouds, fragmented shells and, consequently, new 
generations of clouds. In this work, where we have assumed 
clouds with typical radius rc < hs, the fragmentation time 
due to these interactions is: 



td,sh 



hsh 



[Vsh - Vc) 



(13) 



with Vc and Vsh being the clouds and the shell velocity, 
respectively, which is in agreement with the simulations. 



i:- ■' 





^■l-7- I., 



nf y-. K 



Fig. 4. Color-scale maps showing the time evolution of the 
midplanc density distribution (in log scale) of the inter- 
action of two SNR shells with a system of 4 clouds with 
mc = O.OIM0, rc=0.09 pc, and nc=100 cm"^. The shells 
are injected at times t ~ and t = 4.5 xlO^ yr. Each 
shell has a density nsft=0.09 cm~^ and velocity of 250 km 
s~^, and the ambient medium has n=0.01 cm~^ and T= 
lO*' K. The computational box has dimensions 4 pc x 1 
pc X 1 pc, corresponding to 512 x 128 x 128 grid points 
at the highest grid level. Outflow boundaries have been 
assumed. The times depicted are: t = 1.3 xlO'* yr (loft), 
t = 5.7 xlO'' yr (center), and t = 8.6 xlO'^ yr (right). 
The density is shown in units of g cm~^ (see Melioli, de 
Gouveia Dal Pino & Raga 2004 for details). 

where rc is in pc, Tg = T/10^ K, the factor (j) < 1 allows 
for a reduction in the heat conduction due to the presence 
of magnetic fields, turbulence, and other effects, ^ and 



Kinetic energy is also transmitted by SNRs to the 
clouds, which is lost by cloud-cloud c;ollisions. This bal- 
ance is quite difficult to simulate numerically. In our model 
we simply assume that the cloud velocity stays constant 
at the fiducial value 10 km s~^. 



3.5. Thermal evaporation 

A cold and dense cloud surrounded by a tenuous and hot 
gas of number density n and temperature T may undergo 
evaporation at a mass rate given by (Cowie, McKee and 
Ostriker 1981): 



Mev (Mq yr-1) = < 



-4.34 X 10 
-5.92 X 10 



5/2 



1.3 X 10-«T«' r, 



''^(j)Tg^'^rca~ 
,5/2 



-5/8 



ff C7 < 1 
ff C7 > 1 
if (T < <Tc 



1.54xl07K/ nr, 



1 



(15) 



' c,pcv 



is the saturation parameter. For a < o'er clouds condense 
rather than evaporate (McKee & Cowie 1977; McKee & 
Begelman 1990). These latter authors find = 0.028, 
but Borkowski, Balbus & Fristrom (1990) claim that the 
value should be five times larger: a^i = 0.15. 

Also in this case the mass-loss rate can be obtained as 
a function of the cloud mass m.c whith the help of equation 
(8). In the discussion in the next section, we focus on the 
classical case {a < 1), which represents an upper limit to 
the mass evaporation rate. In this case: 



Mev = -6.94 X 10 



.9^y5/2^_ 



1/7^8/21 



M© yr" 



(16) 



^ We notice that in the simulation of Fig. 4, wc have taken 
clouds with densities nc=100 cm~^; for larger densities (as 
those we have considered in § 5) we must expect a similar 
effect, i.e., an efficient fragmentation of the clouds but without 
significant increase of the ambient gas density. 



the present work, duo to the lack of knowledge of the 
magnetic field configuration in the clouds in a SB environment 
and thus of its truly efi^ect upon the heat transport, we have 
for simplicity adopted <j) = 1. Although this will tend to maxi- 
mize the thermal evaporation rate of the clouds, we find that 
for the typical conditions of the clouds in the SBs the evapo- 
ration will, in general, be dominated by the photoevaporation 
ablation process (see § 4 and Fig. 6 below). 
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4. Time scales 

The energy released by the SNRs and the evolution of the 
ISM will depend on the shells interaction time, Unt- 

A SNR will form only after the SN shock front enter 
the Sedov phase. When the SN explodes, the star ejects 
a mass Mf,j ~ 10 Mq into the ISM with a terminal ve- 
locity ~ 10"^ km s~^ and the ejecta will expand at nearly 
constant velocity until they encounter a comparable mass 
of ambient medium. This occurs at a time t^h which de- 
termines the onset of the SNR formation (e.g., McCray 
1985): 



tsh — 



200 ,M, 



ll/3 ^ M, 



fi Ni/3 



(17) 



where n is ambient number density. The present model 
will be valid only for ts}i < Unt- 

To compute the shells interaction time, tint, we will 
borrow the porosity concept first introduced to investi- 
gate the possibility that the collective effects of the SNRs 
in our own Galaxy might be capable of generating and 
maintaining a rather widespread component of much hot- 
ter, lower density gas in the ISM (Cox & Smith 1974). 
Here, we will assume that tint corresponds to a porosity 
value (Cox & Smith 1974, McKee, 1992): 

TZVsNtint 



Vsi 



I 



(18) 



where Vs n is the volume of a SNR at the time tint and Vs b 
is the volume of the SB ambient. Q gives, approximately, 
a measure of the hot gas filling factor in space-time. When 
it is close to unity, all the ambient medium must be filled 
by SNRs. The shells collide with each other and a unique 
superbubble of hot gas will fill the SB region and may 
drive a galactic superwind. ^ 

Shells interactions can occur before or after the SNRs 
become radiative. Depending on this, we have two possible 
expressions for tint ■ If the shells interact when the SNRs 
are still in the Sedov phase: 

,100 pc.i^/" 



9 X lO'' ( 



x( 



RSB 

3.3 X 10-4 



yr- 



5/11 



n 3/11 



yr (19) 



and if the shells interact when the SNRs are already in 
the radiative phase, then: 



t 



int—rad 



inn nc 21/13 
= 7.5x10' (i^) 



B 



^ Although there has been some debate in the literature re- 
garding the uncertainty in the determination of the value of the 
porosity induced by the remnant population in the ISM of our 
Galaxy and its real importance in providing and maintaining 
a hot gas component (Cox 1995, McKee 1995; see also Slavin 
& Cox 1992, 1993, McKee & Zweibel 1995), in the case of SB 
galaxies, the SN rate is more than 10 times higher than in nor- 
mal galaxies and therefore, SNRs interactions are expected to 
be efficient enough to eventually provide a substantial value 
for the porosity. 



, 1 1 1 1 1 1 1 


, , 1 1 1 1 1 1 1 1 1 


t 




'mt-rad 













log n 

Fig. 5. Time scales versus gas ambient density of a SNR in 
a spherical SB region, with Rsb = WO pc. The interaction 
time in the Sedov phase is represented by the dotted line, 
the interaction time in the radiative phase by the dashed 
line, the time scale for the radiative cooling by the solid 
line, and the time scale for the formation of the shell by 
the dot-dashed line. 



3.3x10-4 yr-\7/i3 



/3° 



£^51 



(20) 



where P has been defined in § 2.2 . Figure 5 compares 
these interaction times (as functions of the ambient den- 
sity) with the time scale for the shell formation, tsh, and 
with the radiative cooling time of the SNR shell (McCray 
1985): 



t, = 1.49 X 104/3-5/14^-4/7 



(21) 



We notice that for ambient densities < n.l cm ^, we must 
take tint = tint-sedov, siuce tint-sedov < tc m this range. 
For t > tc, the SNR enters the radiative phase and tint-rad 
should be employed. Since for most of the calculations 
of interest in this study (see below), we take either n = 
0.01 cm-3 or 5 cm-^, we will use tint = tint-sedov or 
tint = tint-rad, respectively, to evaluate the time scale for 
superbubble formation. 

Using equations 7, 12 and 16, and defining the time 
scale for the ambient density increase as the time a given 
cloud mass-loss process (photoevaporation, thermal evap- 
oration or drag) takes to increase the ambient number 
density by 1 cm-^, then we have: 



tuv 



ggy^ RsB ^3(8 X 106 MQ,me3/7 



100 pc' 



Ma 



tev = 1.9 X 10'^( 



Rsi 



100 pc 



n- 



X 10^ 



yr 



rji 5/2 



(22) 



yr(23) 
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2/3 1/3 



yr(24) 



for the density increase timescale through photoevapora- 
tion, thermal evaporation and classical drag, respectively. 
In the equations above, n and ric are given in cm~^. 
Depending on how all these time scales relate to each 
other, one finds either a superbubble scenario or a gas 
mixing scenario. When tint is longer than one of the three 
time scales above, the mixing of the ablated matter from 
the clouds occurs before the formation of a high pressure 
bubble. In this case, the ablated gas will enhance the den- 
sity of the hot diffuse gas of the ISM and increase its 
radiative cooling. In contrast, if Unt is shorter than the 
destruction time scales above, a low density superbubble 
will form very fast and the SN heating efhciency (HE) will 
rapidly increase to unity. 

Figures 6 compares these time scales for typical param- 
eters of a SB region, where the time scales are depicted as 
a function of the mass of the clouds. We see that for the 
two adopted pair of values of the density and temperature 
of the ambient medium, n = 0.01 cm~^ and T = 2 x 10^ 
K, and n = 5 cm~'^ and temperature T ~ 5 x 10'* K, the 
presence of clouds with mass TOc < 100 Mq will make tint 
longer than at least one of the time scales for clouds abla- 
tion. In this case, the gas mixing scenario will be dominant 
and the formation of a superbubble is postponed. 




5. An evolution model for HE 

5.1. Equations 

Once considered all the relevant physical processes, we 
can now write the equations governing the evolution of 
our model described in section 2.1. In particular, we de- 
scribe the time evolution of the ambient gas density, n, 
and pressure, p, and of the clouds number Nc with a mass 
TUc. The set of equations is given by: 



dMio 



dt 
dt 



Muv + Mev + Md (Mq yr-1) 



Ml, 



Ci- TVc 



-C2N,v,-C^N,-^ (yr"')(26) 

tcoll 



dfl • TL 

— = CaNcMioss - C5- C2nMcs 




m, (M, 



(25) Fig. 6. Steady time scales versus mass of the clouds in a 
spherical SB region, with i?sB=100 pc, J^=0.225, Mg = 
8 X IO^Mq, -0 = 1. The interaction time is represented 
by the dot-dashed lines, the ambient density increase 
timescale due to photoevaporation by solid lines, due to 
thermal evaporation by dotted lines, and due to drag by 
dashed lines. Upper panel: n=O.Gl cm"'^, T=2 x 10^ K; 
bottom panel: n=5 cm^"^, T=5 x 10* K. 



+ C6(^) (cm-^s-i) 



CiopMcs (erg s ^) 



(27) 



(28) 



Eq. (25) computes the total mass-loss rate per cloud us- 
ing the equations of the clouds mass-loss rate (7), (12), 
and (16); Eq. (26) gives the evolution of the total number 



of clouds; and Eqs. (27) and (28) give the mass and en- 
ergy conservation, respectively, of the diffuse gas. In these 
equations, Cn (with n = 1 to 10) are constants which are 
given in Appendix B, Cg is the sound speed of the ambient 
gas, A4 its Mach number, that is taken to be equal to 1 
in this study, and Vc is the velocity of the clouds which is 
assumed to be 10^ cm s^*, as discussed in § 3.4. In Eq. 
(26), the first two terms give the net difference between the 
processes of cloud formation and cloud ablation, the third 
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term gives the rate of clouds escaping through the bound- 
aries, and in the last term we have introduced the typical 
time scale for cloud-cloud collisions, tcou = l/(-^c47rr^Wc), 
in order to account for the possibility of collisions between 
clouds in the SB volume. We note that we have assumed 
a constant mass for the clouds, rric. This is equivalent to 
adopting a constant mass distribution function (or a con- 
stant average mass), and is a valid assumption, as long as 
the processes for cloud formation and destruction do not 
affect much the average mass of the clouds but only their 
total number. 

In Eq. (27), the first term gives the ambient gas den- 
sity increase due to the total mass-loss rate by the clouds; 
the second term gives the density decrease due to the in- 
teraction with SNR shells that sweep the ambient mate- 
rial; the third term represents the gas escaping through 
the boundaries, and Ce gives the injection of matter from 
stellar winds and SNe explosions (Leitherer & Heckman 
1995). 

Eq. (28) computes the energy of the diffuse gas. The 
first term is the rate of energy injection from the SNe, 
the second gives the kinetic energy term of the gas lost 
by the clouds, the third gives the energy decrease due to 
the radiative cooling and the last term gives the energy 
density decrease due to the escape of the gas through the 
boundaries. 

Using a fourth-order Runge-Kutta integrator, we de- 
termine the evolution of the SN efficiency HE defined as 
(see section 2.1): 



HE 



47ri?2 



(1 



pv 



SB' 



TIE 



SN 



(29) 



5.2. Results 

In sections 2 and 3, we have studied the physical processes 
that drive the mass- loss rate of the clouds and these have 
been compared with the process for superbubble forma- 
tion. We have found that, while the time scale for the in- 
crease of the ambient gas density, tmx — rnin{tijv Ttevitd), 
is shorter than the time scale for the superbubble expan- 
sion, tint^ the HE value must be small. When, on the other 
hand, tint < tmx, then HE must become close to unity. If 
this is correct, one should expect an increase in the HE 
value with the decrease of the ambient gas radiative cool- 
ing, which is oc n^. 

In Figure 7, we see the results for a reference model 
that we have called Model 1 . We have here considered neu- 
tral clouds (Tc^lOO K) which are initially in total pressure 
equilibrium with both the ambient gas and the ionizing 
photons. Table 2 summarizes the input parameters for this 
model. The SB has initial conditions which are taken from 
the limits given in Table 1. We notice that the SN heating 
efficiency increases to 1 after 16 Myrs. Before this time, 
the total energy stored in the gas is only a few percent of 
the total energy released by the SNe. Most of it is radi- 
ated away, as a result of the heating of the dense gas. The 
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Fig. 7. Model 1: Ambient density (n), temperature (T) 
and pressure (p) evolution in the top-left panel for a spher- 
ical SB with initial Mb = 10^ M©, i?sB = 100 pc, m = 0.01 
cm-3, To = 2 X 10^ K, po = 6.35 x lO^^^ cm^^ and 
total gas mass Mg = 8 x lO^M©. In the bottom- left panel: 
number of clouds (Nc) and radius of clouds (rc) evolution; 
in the top-right panel: the SN heating efficiency (HE); 
and in the bottom-right panel: the clouds filling factor, 
(// — NcVc/VsB, where Vc is the volume of a single 
cloud), the saturation parameter (cr, eq. 15) and the cloud 
to ambient gas density ratio (x — Pc/p). The ambient 
density, temperature and pressure, and the number and 
radius of clouds are normalized to their initial values (in 
log scale). The initial radius, temperature and density as- 
sumed for the clouds are Tc = 6.5 x 10"'^ pc, Tc = 100 K, 
and ric = 6 X 10^ cm"^, respectively. 



mass-loss rate from the clouds due to the processes of pho- 
toevaporation, drag and thermal evaporation for Model 1 
is shown in Figure 8. The domination of the photoevap- 
oration during the first 16 Myrs is clear, and highlights 
the fact that after 16 Myr no clouds are present in the SB 
(see bottom- left panel of Figure 7). Thus, after this time, 
the clouds mass-loss becomes unimportant for the physi- 
cal evolution of the system. We also note that in this case, 
instead of thermal evaporation the clouds actually suffer 
a thermal condensation, through M^y (as cr < Cc = 0.03, 
in the bottom-right panel of Figure 7), but this does not 
influence the global results of the HE evolution. 

The initial large increase of the density (by four or- 
ders of magnitude) in Figure 7, may seem unrealistic, in 
principle, but it is actually due to the choice of the ini- 
tial conditions. Since we have run our model with an ISM 
with initial conditions which are very similar to those of 
a superbubble (n = 0.01 cm^^^ T = 2 x 10^ K) and 
with tmx < tint, Figure 7 shows that the ISM density in- 
creases almost instantaneously to a typical value for SBs 
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Table 2. Values for the parameters of the models assumed in this work 
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i (Kj 




Cloud mass-loss 






1 


0.01 


2 X 10'' 


0.015 


MuV, Mev, Md 


8 X 10^ 




2 


5 


5 X 10" 


0.015 


Muv, Mcv, Md 


8 X 10*^ 


10*^ 


3 


0.01 


2 X lO'* 


0.015 


MuV, Mev, Md 


8 X 10^ 


10^ 


4 


0.01 


2 X lO'* 


0.015 


Mev, Md 


8 X 10*^ 


10® 


5 


0.01 


2 X 10^ 


mass dist. func. < m > = 0.3 


Muv, Mev, Md 


8 X 10® 


10® 


6 


0.01 


2 X lO'' 


0.015, no clouds formation 


Muv, Mev, Md 


8 X lO'' 


10® 




Fig. 8. Model 1. Mass-loss rate of the clouds with mass 
TOc = 0.015 M0 due to photoionization, Muv (solid hue), 
thermal evaporation, Me„ (dotted hue), and drag, Md 
(dashed line). 



Fig. 9. Model 2. The same as in Model 1 of Figure 7, 
except that the initial ambient gas density and tem- 
perature are rt=5 cm^'^, T = 5 x lO"' K, respectively 
(po = 3.45 X 10-"). 



(n ~ 10 — 50 cm"-^). Therefore, the assumed initial su- 
perbubble disappears before being able to blow out the 
gas. Indeed, if we run a model with initial conditions that 
are typical of a SB, i.e., with ambient T = 5 x 10^ K 
and n = 5 cm~^ (see Figure 9, Model 2 of Table 2), the 
density increases only by an order of magnitude, but the 
corresponding HE value remains small during the same 
period of time. 

The comparison of the results of Figure 7 and 9 above 
indicates that the HE evolution is very similar in both 
models, or, in other words, it is very insensitive to the 
initial conditions of the hot ISM. 

In both, as long as the ambient gas density remains 
large enough for the radiative cooling to be efficient, the 
HE value remains very small. After 16 Myr, when 
the density of the gas available in the SB becomes too 
small (< 10~^cm-'^) due to continuous escape through 
the boundaries of the system and exhaustion of the gas in 
the clouds, then the radiative cooling ceases and a sudden 
increase in the temperature leads to a superbubble forma- 



tion in a time scale Unt — 10^ yrs. In consequence, HE 
1. 

In Figure 10, we have run a model (Model 3 of Table 2) 
with the same parameters of the Model 1 , but with a total 
stellar mass 10 times bigger. In our model, this implies an 
increase in the SN rate explosion (7^) and the total gas 
mass in the clouds [Mg) of a factor 10. Under these con- 
ditions, we see that HE never reaches high values during 
the SB lifetime due to the larger gas mass in the clouds. 
After an initial increase caused by the first explosion of 
the SNe at a higher rate, the HE value decays to less than 
0.5 per cent. 

In Figure 11, Model 4 presents results for the HE evo- 
lution considering clouds which are fully ionized (Tc = 10^ 
K). In this case, the photoevaporation was not taken into 
account, since the clouds are transparent to the ionizing 
photons. As a consequence, in this case they are in pres- 
sure equilibrium only with the ambient medium. The com- 
parison with Model 1 of Figure 7 shows that the HE value 
increases slower in Figure 11 (where HE ^ 1 at 20 Myrs), 
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because in the case of Figure 7 the clouds are more ef- 
ficiently destroyed, and survive for a smaller time in the 
SB. 

At this point, it is interesting to highlight that the 
clouds are the real valves that regulate the feeding of the 
diffuse ambient medium of the SB, in the sense that they 
are able to simultaneously retain part of the gas to them- 
selves and to lose part of it to the ISM through cloud 
destruction, thus maintaining the ISM density at high val- 
ues, as long as tmx < tint- 

Using a cloud number distribution (which is given in 
Appendix A) N{mc) oc m-^", with a = 1.5 and an av- 
erage cloud mass < m >~ 0.3 Mq (Model 5) instead of 
a single initial cloud mass, the previous results of mod- 
els 1, 2 and 4 (for which a single cloud mass — 0.015 
Mq was adopted) are not very much changed, as we see 
in Figure 12. In this case, the density grows by 3 orders 
of magnitude at the beginning, and the number of clouds 
evolve more smoothly than in the other models due to the 
presence of a cloud mass distribution instead of a single 
mass. Nonetheless, since the density evolution is driven 
mainly by the total amount of gas in the clouds, which is 
the same as in the other models, the HE value goes to 1 
after ~ 16-19 Myr, i.e., just like in the previous models. 
This result also shows that HE is quite insensitive to the 
adopted cloud mass, rric, since the average mass of Model 
5 is almost 10 times larger than the adopted cloud mass in 
the previous models {rric — 0.015 Mq). Consistently, tests 
made with a constant cloud mass tUc equal to the average 
mass of the distribution function above (< rric >~ 0.3 
Mq), have produced similar results to those of Model 1 of 
Figure 7. 

All the models investigated above (Figs. 7 to 12), have 
assumed a continuous formation of new clouds, during 
the SB lifetime, which are generated by swept ambient 
gas accumulated on the SNR shells. If we now consider 
an extreme situation in which no deposited gas on the 
shells is allowed to form clouds (see model 6 of Table 2), 
then in this case Figure 13 shows that HE quickly goes 
to one in about 2 Myr. This situation is equivalent to as- 
sume a fully ionized gas both in the shells and the ambient 
medium with negligible radiative cooling, so that the time 
scale for clouds formation from fragmentation of cooled 
material deposited on the shells is much longer than the 
time scale for SNR interactions. In Figure 13 (model 6), 
only the first initial population of cold clouds is present 
in the SB environment and it is soon ablated (mainly by 
photoevaporation) adding cold, dense gas to the ambient 
medium. However, this is not sufficient to prevent the fast 
overlap of the SNRs to form a hot superbubble and the 
temperature increase of the diffuse ambient gas to about 
10^ K. Figure 13 actually gives a lower limit for the time 
it would take for HE to increase to one. In a more real- 
istic scenario, at least part of the deposited material on 
the shells should form new generations of clouds, and HE 
would then increase to unity only after several Myrs, like 
in Figs. 7 to 12. We note further that, if we make all the 
gas, including that of the clouds, initially fully ionized (like 
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Fig. 10. Model 3. The same as in Model 1 of Figure 7, 
except that Mb = lO"^ Mq and Mg = 8 x lO"^ Mq. 
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Fig. 11. Model 4. The same as in Model 1 of Figure 7, ex- 
cept that the clouds are fully ionized (Tc — 10^ K, rc=1.7 
pc, and nc=2 cm~'^), and photoevaporation has been ne- 
glected. 

in Figure 11), so that no photoevaporation occurs after- 
wards and do not allow for formation of new clouds, like in 
Figure 13, then in this case, we obtain that the HE value 
goes to one only after ~ 8 Myr. 

In summary, the results above show that, provided that 
a continuous formation of clouds due to deposited material 
on the SNR shells occurs in the SB, then the common 
assumption of a HE value close to 1 is reasonable only 
after many Myr, when a large quantity of matter has been 
blown-out and the ambient gas density has dropped to 
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Fig. 13. Model 6. The same as in Figure 7 (Model 1), 
except that only the first generation of clouds is present in 
the SB ambient medium, and no further clouds are allowed 
to form. 

values of < 10~^ cm~^ (Note that in Figs. 7 to 13, the 
ambient densities n have been normalized to their initial 
values which are given in Table 2.). 

6. Discussion and Conclusions 

The high stellar density in a SB galaxy allows for SNRs 
interactions that very quickly may fill all the ambient. 



forming a superbubble containing hot gas with tempera- 
tures T=10''-10S K, and densities n=0. 1-0.01 cm'^. Shells 
interactions increase the process of clouds formation with 
typical temperatures in the range of 100-10"* K, depending 
on whether they are photoionizated or not. Under these 
conditions, hot gas may eventually escape from the galaxy 
in an outflow that changes the chemical and dynamical 
characteristics of the galaxy itself. It is clear that single 
SN explosions have a typical efficiency in the energization 
of the ISM of normal galaxies of few per cent, but a de- 
bate exists about the SN heating efficiency in SB galaxies. 
To solve this problem, we have developed a simple semi- 
analytical model that analyses the most important physi- 
cal processes in a three-component system formed by hot 
gas, SNRs and cold (or warm) clouds. Under this condi- 
tion, the evolution of the three-phase medium after the 
development of an instantaneous SB has been examined 
to compute the SN heating efficiency (HE). The study has 
also been accompanied by 3-D numerical experiments in- 
volving the interaction between SNRs and clouds in a SB 
environment. We have found (see § 5) that HE remains 
very low as long as the mass-loss time scale of the clouds 
{t7nx) remains smaller than the time scale for a superbub- 
ble formation (tint)- For an environment with continuous 
formation of clouds, this occurs during the first 16 Myrs 
for an initial total gas mass in the clouds of 8 xIO^Mq, 
after which HE rapidly increases to unity. This result can 
be explained by the fact that the efficient ablation of the 
clouds during the first 16 Myrs increases the density and 
the radiative cooling of the gas of the ambient medium, 
therefore preventing a temperature enhancement (HE ~ 
1) that would lead to gas expansion with the formation of 
a hot and rarefied superbubble and a sudden expulsion of 
the gas through a galactic wind. 

This result is quite insensitive to the initial conditions 
of the hot diffuse ambient ISM of the SB, nonetheless the 
semi-analytical calculations show that the SN heating effi- 
ciency can be affected by parameters such as the total gas 
mass in the clouds. We find that an increase in this pa- 
rameter results low HE values for longer times. In fact, an 
increase in this parameter by a factor 10 with respect to 
the value above makes HE < 1 over the entire SB lifetime. 

We have also found that under an extreme condition 
in which no further clouds form after the first generation 
of clouds, then this initial population alone is unable to 
prevent a faster raise of the ambient gas temperature and 
therefore, of , the HE value to one. However, such situation 
could be expected to occur only in environments which are 
initially already completely ionized by the UV photons 
from massive stars, having no cold clouds or shells. 

In summary, we can conclude that the SN efficiency 
has a time-dependent trend which is sensitive to the ini- 
tial conditions of the SB system. Therefore, it is not pos- 
sible to assume a single value for HE over the entire SB 
lifetime without considering the physical interactions that 
can change the environmental characteristics. If all the 
gas that is swept by the SNR shells goes into clouds for- 
mation, providing a continuous generation of clouds, then 
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HE remains small up 16 Mrys. If, on the other hand, the 
material that is swept by the SNRs does not lead to new 
generations of clouds, then the gas forms a supcrbubbic 
in few Myrs. 16 Myr is, therefore, an upper limit for the 
time it takes for HE to increase to unity. 

Our model allows one to obtain more realistic HE val- 
ues depending on the evolution stage of a SB galaxy. In 
particular, it provides a natural explanation for the low 
values invoked in some models to describe the chemical 
evolution of dwarf galaxies (Carigi et al. 1995; Bradamante 
Matteucci & D'Ercole, 1998; This is the case, for example, 
of the SB galaxy SBIZwlS (Recchi, Matteucci & D'Ercole 
2001), for which a value of HE 0.03 is required to ex- 
plain observations and is in good agreement with the re- 
sults of this work which predicts HE ^ 0.01 to 0.1 in the 
first < 16 Myr of the SB lifetime (Figs. 7-12). 

We have conducted our study considering regions 
which are small compared to the total size of the host 
galaxy. There must be other relevant physical phenomena 
outside of this region that can also influence its evolution 
and the determination of HE, as for example, the rate at 
which material is accreted from the neighborhood into the 
SB region. In fact, during some period of time, SBs may 
be effectively fed by gas from outside. The 2D hydrody- 
namic models of galactic winds in SB galaxies found in 
the literature often produce a superbubble that exceeds 
the dimensions of the narrow cones detected in the nu- 
clear SB, but Tenorio-Tagle & Mufioz-Tufion (1998) have 
shown that, if the gas of the galactic disk is accreted onto 
the SB region at a rate of a few Mq yr~^, the funnels 
may present small opening angles, in agreement with the 
observations. This infall is thought to be necessary to ac- 
count for the accumulation of the matter detected in these 
sources (Sanders & Mirabel 1996) and the large rates of 
star formation (Larson 1987 et al., Suchkov et al. 1994). 
In the present study, although these effects of matter ac- 
cretion from the neighborhood have not been taken into 
account, we expect that their presence will keep the HE 
value smaller than unity for even larger times, since they 
must cause an increase in the formation rate of clouds. 

Recently, IR, radio, continuum and UV observations of 
SBs have revealed the presence of very compact structures 
with typical sizes of few parsecs and mass in the range of 
10^ — 1O^M0, that can be associated with star clusters 
and super star clusters (see, e.g., Ho 1997; Johnson et al. 
2001; Gorjian, Turner & Back 2001). This is somewhat 
different from the general assumption often employed in 
the modeling of SBs and also adopted in the present work, 
that the star formation rate is diluted over all the SB vol- 
ume. However, even if the bursts are not homogeneously 
distributed, we expect that the results here obtained will 
not be much affected. Since each star cluster must drive 
a bubble with morphological characteristics similar to a 
SNR, when these bubbles interact with each other the 
global evolution within the SB will be similar to the one 
analyzed in the present work. 

Another relevant aspect that deserves attention is the 
determination of the conditions of the galactic wind that 



develops when HE ~ 1. In this work we have assumed a 
constant rate for the SN energy injection over the SB life- 
time. This assumption is consistent with previous mod- 
eling of SB systems (see, e.g., Leitherer et al. 1999). 
However, if we had taken into account another possible 
evolution scenario, with a decreasing SN energy injection 
rate with time, then in this case, the reduced available SN 
injection power in the late stages could be not suflacient 
to drive the galactic wind that, according to our model, 
should evolve only in the second half of the SB lifetime. 
This may offer a potential explanation of why some SB 
galaxies produce winds and others do not. The computa- 
tion of the total quantity of mass that is transported by 
the galactic wind (Suchkov et al. 1994) is also dependent 
of the physical history of the expelled clouds. Although at 
a different rate, they will be also ablated by photoevapo- 
ration, drag and thermal evaporation, contributing to the 
growth of the matter that forms the outflow itself. 

A more particular attention should be given to mag- 
netic fields, which can be important in reducing the abla- 
tion processes and which can influence the escape velocity 
of the hot gas and the cold clouds. Although in the present 
work they have been taken into account in the computa- 
tion of the total internal pressure of the clouds (eq. 6) and 
also explicitly in the evaluation of the their dominant ab- 
lation process (the photoevaporation), they have not been 
included in the evaluation of the pressure of the ambient 
gas and the SNRs shells, or in the numerical simulations of 
SNR-clouds interactions. Also, to be remarked is the fact 
that the adopted value, of ~ few /iG within the clouds, is 
somewhat uncertain since, to our knowledge, there is no 
direct magnetic fleld measurements in SBs. Nonetheless, 
it is consistent with the values often required for pressure 
equilibrium of the clouds with the hot ambient medium in 
our Galaxy, and also with the lower limit strengths nec- 
essary to collimate the galactic winds in SBs (de Gouvcia 
Dal Pino & Medina Tanco 1999). In the case of the SNRs, 
while the presence of magnetic fields in the shells could, 
on one side, reduce their compressibility and thus decrease 
the volume of hot gas within the remnant (see. e.g., Slavin 
& Cox 1993, McKee & Zwcibcl 1995), on the other side, it 
could constrain the thermal conduction and thus inhibit 
the destruction of this hot gas component (e.g., McKee 
1995). The reduction of compression by the magnetic field 
could also inhibit the process of shell fragmentation dur- 
ing interactions with other SNRs and, consequently, the 
formation of clouds by this process (that was discussed in 
§ 3.1). The relative importance of these competing effects 
may be quantified after determination, from observations, 
of the magnitude and geometry of the magnetic fields in 
SB environments. 

Finally, we should notice that studies that search for a 
potential connection between SBs and AGN activity, that 
suggest that SBs at the nuclear regions of active galaxies 
could eventually trigger the development of a supermas- 
sive black hole (BH) in the center, may be favored by the 
results of the present analysis as these predict a delay in 
the expulsion of gas from the SB in form of a galactic 
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Fig. A.l. Initial number of clouds for a distribution of the 
clouds number in a range of 10-"^ - IO^Mq, and the value 
of the average mass < m > of the distribution Nc oc m~", 
with a = 1.5 (left), and a ~ 1.9 (right). 

wind, for ~ 16 Myrs. The gas retained in the system may 
lead to new bursts of star formation and increasing depo- 
sition of material in the center of the SB that, in turn, 
may eventually form a BH. 
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Appendix A: 

The cumulative cloud mass-loss rate in a SB may be com- 
puted taking into account the spectrum of the clouds 
number N{m) cc m~" between an upper (toc,sup) and 
a lower (mc,inf) limit, instead of considering only a single 
cloud mass. Observations of high star formation regions 
indicate a minimum mass of 10~*-Mq and a maximum 
mass of ^ IQ^Mq (Kramer et al. 1998). This is consis- 
tent with the theoretical maximum stable mass of a mag- 
netically dominated cloud against gravitational collapse 
(Bertoldi & McKee 1990). We thus assume TOc,sup = 10^ 
Mq and X = mc.inf/^c.sup = 10^^. Also, from the stud- 
ies of Zinnecker, McCaughrean & Wilking (1993) and Blitz 
(1993), we take a — 1.5. A plot for this cloud mass distri- 
bution is shown in Figure Al. 

Let /{rric) = 6m2 be the mass- loss-rate of a single 
cloud, where the constants 9 and 7 depend on the spe- 
cific process responsible for the cloud mass-loss studied 
in § 3 (i.e. thermal evaporation, drag or photoevapora- 
tion). The cumulative mass- loss rate must be computed 




Fig. A. 2. The function g(a,j,X) for a — 1.5 (left panel) 
and a = 1.9 (right panel) for three different values of 7. 



taking into account the mass spectrum of the clouds above 
dNc{mc) / dm oc m^T". The total mass- loss rate per process 
is thus given by 

/■nicsup 

M = / Nc{mc)f{mc) dmc, 

J mc,inf 

and can be written as 

M ^ NJ{m,^snp)gia,-f,X) 

where Nc is the total number of clouds inside the SB region 
and 

5 a,7,A = — . 

1+7 — a 1 — " 

Taking into account Equations (7), (12) and (16) for Mjjv, 
Md and Mev respectively, we obtain 7 ~ 12/21, 16/21 and 
8/21 for cloud photoevaporation, drag and thermal evap- 
oration processes, respectively. We find that g{a,"f,X) 
is rather sensitive to the parameters involved. Figure 14 
shows 17(0;, 7, X) as a function of A", for a = 1.5, a = 1.9, 
and the three values of 7. 

For T = 0.225 (see eq. 6), we obtain Muv = 2.9 x 
lO'^Mg Mq yr-i, Mcv = 6.17 x IQ-^^T^'^Mg M© yv'^ , 
and Md = 1.32 x lO-^^nx^^^we^g Mq yr~^ As demon- 
strated in Section 3, photoevaporation is by far the most 
important mechanism for the cloud mass-loss rate during 
the first Myrs of the SB life, even considering a cloud mass 
distribution, and the results do not differ much from those 
when one takes a single initial mass (see section 5). 

Appendix B: 

We present in this Appendix the explicit expressions for 
the constants Cn, with n = 1 to 10, that appear in equa- 
tions 25, 26, 27, 28 in § 5. 

4.64 X 10^4 FsamHM ,^ 
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where S is the fraction of gas that is accumulated in the 
SNR sheU. 



a 



1 



3.15 X 107 rur 



4.6 X 1029 VSB 
where Vss is the SB volume 



C5 = 



VsB (3.08 X 1018) 



(B.3) 



(B.4) 



(B.5) 



where Ssb is the spherical external surface of the SB in 
pc. 



1.1 X 10-24 

VsBmafJ- 



(B.6) 



which gives the increase in the density due to the mass-loss 
rate from the stellar winds and the SNe explosions. 



C7 = 



Co 



c 



10 



7.24 X 10-6* 
4.56 X 10-23 



'2SsB 

VsB 



(B.7) 
(B.8) 
(B.9) 
(B.IO) 



All the dimensional quantities of the SB 
{VsB, Ssb, Rsb) are expressed in pc. The mass-loss 
rate of the clouds is in Mq/jy, the time scales are in yrs, 

0- 



and we rric, that is assumed constant, in MfT,. 
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